Theory of strong-field injection and control of photocurrent in dielectrics and 

wide bandgap semiconductors 
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We propose a theory of optically-induced currents in dielectrics and wide-gap semiconductors 
exposed to a non-resonant ultrashort laser pulse with a stabilized carrier-envelope (CE) phase. In 
order to describe strong-field electron dynamics, equations for density matrix have been solved self- 
consistently with equations for the macroscopic electric field inside the medium, which we model 
by a one-dimensional potential. We provide a detailed analysis of physically important quantities 
(band populations, macroscopic polarization, and transferred charge), which reveals that carrier- 
envelope phase control of the electric current can be interpreted as a result of quantum-mechanical 
interference of multiphoton excitation channels. Our numerical results are in good agreement with 
experimental data. 
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I. INTRODUCTION 



The time it takes electrons in a solid to respond to 
an external electric field is on the order of a fraction 
of femtosecond. Due to recent progress in ultrafast 
and attosecond science, it is nowadays possible to per- 
form attosecond-scale time-resolved measurements of ex- 
tremely nonlinear phenomena that occur when intense 
few-cycle laser pulses interact with a solid, and oppor- 
tunities related to this progress have recently begun to 
be explored. In particular, implicit observations of sub- 
cycle temporal structures associated with interband tun- 
neling in a dielectric were reported^'^, and effects related 
to Bloch oscillations in a bulk solid were observed^"^. 

A direct observation of attosecond-scale electron mo- 
tion was reported by Schiffrin et al. in a very recent pa- 
per^ demonstrating that electric currents in a fused silica 
sample can be switched and driven by the instantaneous 
field of an optical waveform at intensities just below the 
damage threshold. It was argued that, at field strengths 
where the induced potential difference between neighbor- 
ing unit cells approaches the bandgap energy, the exper- 
iment can be interpreted in the basis of Wannier-Stark 
states. In this paper, we show that this kind of measure- 
ments can also be explained within a more conventional 
approach based on interference between different multi- 
photon channels. This relates the results of Ref. to the 
scope of coherent control. 

Injection and coherent control of electric currents in 
semiconductors was studied before, both in experiments 
and theory^"— . Coherent control was demonstrated in (w, 
2a;) schemes based on the interference of one- and two- 
photon excitation pathways induced by a laser pulse and 
its second harmonic. In such experiments, it is sufficient 
for pulse durations and time delays to be less or compa- 
rable to the carrier dephasing time 100 fs). Now, with 
the availability of the few-cycle pulses, control over the 



photocurrent can be achieved within a single laser pulse 
and on much shorter timescales ('^ 1 fs). 

The interpretation of carrier-envelope phase effects in 
atomic and molecular systems in terms of quantum- 
mechanical interference is also well know n^'^i^^ , and a 
general abstract theory that describes all phase effects 
in these terms is availablei^. However, it is not obvi- 
ous whether concepts developed for relatively weak fields 
can be applied to interpret experiments with intense few- 
cycle pulses, in which the absorption of many photons is 
required to excite valence-band electrons, and laser field 
amplitude reaches values at which perturbation theory is 
expected to break down. 

To investigate the non-resonant optical injection and 
control of electronic currents in dielectrics, we develop 
a model based on a self-consistent solution of multiband 
optical Bloch equations (OBE) ^'^'^^ together with equa- 
tions for the dielectric polarization and field inside the 
crystal. The paper is organized as follows. In Sec. HIl we 
discuss our theoretical formalism, the gauge choice, and 
approximations that we make calculating the amount of 
charge flowed through a capacitor-like junction. Results 
of our numerical simulations, their interpretation and 
comparison to experimental data are given in Sec. IIIII 
Sec. IIVI presents our conclusions. 



II. THEORY 

In order to calculate the current and polarization in- 
duced by an ultrashort pulse, we solve the density- 
matrix equations in the independent particle approxi- 
mation. Unlike the time-dependent Schrodinger equa- 
tioriF^ii^, this approach takes into account Pauli block- 
ing of interband transitionsil and it can be extended 
to account for electron-electron scattering and interac- 
tion with the bath. The approximation of independent 
particles is reasonable for ultrafast strong-field phenom- 
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ena since the carrier-field interaction is much stronger 
than the carrier-carrier interaction. Another argument 
in favor of this approximation is that the net current 
only depends on the total electronic momentum, which 
is not changed by electron-electron interactior^i^iiS. On 
timescales much shorter than a period of lattice oscil- 
lations, we can also neglect the electron-phonon inter- 
action. For LO-phonons in typical semiconductors and 
dielectrics, the oscillation period is about tens of fem- 
toseconds. 

In the basis of Bloch states, the Hamiltonian of the 
electronic subsystem of a dielectric interacting with a 
laser field can be written as 



H = Hn 
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where £ = {/, k', s'} and j — {i, — k, s} are composite 
indices that include a band index {i and /), crystal mo- 
mentum (k and k') and spin (s and s'). 

In the velocity gauge, the amplitude of interaction with 
the optical field can be written as 



TWyf (k',k,i) = -a/,— A(t)^(k'-k)[%fik + p/,(k)] 
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where A(i) is the vector potential in a medium, indices 
i and / enumerate all valence and conduction bands, Ufi 
is a band-specific sign function defined as 
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is the momentum matrix element between Bloch ampli- 
tudes Wi^k, ^ being the volume of an elementary cell. 

The term proportional to A'^(t), can be eliminated in 
the dipole approximation by the following gauge trans- 
formation^S: 



<if'{t) = *(<)exp 
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The expression for the interaction amplitude in the 
length gauge can be written as^liSS 

Mhf{k',-k,t) = -a/,eE(t)[%iVk + 0^(k)]'5(k' - k). 



where E(t) is the electric field in the medium, and 



= ^ / dVu}_kWVk«.k(r). 
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In the present paper, we use the velocity gauge be- 
cause it results in a system of dynamic equations where 
different k-states are uncoupled. In contrast, the cor- 
responding length-gauge equations have a term propor- 
tional to Vk(5(k' — k), which couples different k-points 
and introduces singularities, and matrix elements ^/i(k), 
which are not uniquely defined. These problems with the 
length gauge in the description of optical phenomena in 
solids are well-known, and have been solved in the recent 
theoretical treatment a^^i^^ . 

Nevertheless, the treatment of the crystal polarization 
in the velocity gauge has a few important drawbacks, 
discussed, for example, in Ref. [2^ First, the solution 
of dynamic equations requires a large number of bands, 
only a few of which contain a significant carrier popula- 
tion at the end of a simulation. Second, a certain sum 
rule for momentum matrix elements must be satisfied^^. 
Violation of these conditions leads to serious numerical 
artifacts and unphysical results. In other words, a de- 
scription of strong-field phenomena in the velocity gauge 
requires an accurate solution of the stationary problem. 
Equations in the length gauge have an advantage of less 
strict requirements to the quality of the eigenproblem so- 
lution, so that the energy spectrum and optical matrix 
elements could be calculated with more simple and rough 
methods, or even used as adjustable parameters^". 

Despite these problems, a fully ab initio approach 
based on time-dependent density functional theory in 
the velocity gauge was developed^iiS^ and successfully 
applied to describe the non-perturbative polarization re- 
sponse of solids to a strong field. Such calculations are 
very resource-demanding and require a high-performance 
supercomputer. 

Our present approach is based on the numerical solu- 
tion of density-matrix equations for a one-dimensional 
lattice and has very moderate computational require- 
ments. Yet it allows to reproduce important features of 
laser-matter interaction on the ultrashort timescale for 
intensities near the damage threshold, and it gives good 
agreement with experimental data. 

Starting from the definitions of two-point density ma- 
trix elements 

(k) _ /J \ (k) _ ; t \ 

Pji = {dj^-kCiM), Pee' = (Cf,kC£',k), 
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Eqs. and the Liouville-von Neumann equation 
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one obtains the following system of equations for each k 
in the Brillouin zone: 
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where 



e 



£u' — Swee + Aiui, 



and indices i and j enumerate conduction and valence 
bands, respectively. Quasimomentum indices are omitted 
here, for simplicity. 

The charge carriers generated by an intense laser pulse 
induce a strongly nonlinear response. If the characteristic 
size of the sample in the direction of polarization is small 
enough, then the field inside the crystal is affected by the 
surface charge distribution, which should be taken into 
account by solution of semiconductor Maxwell-Bloch 
equationsHiSSiiSI. However, a full numerical solution of 
these equations in a three-dimensional case presents a 
rather hard computational problem, so we resort to a 
simplified model of the dielectric polarization. 

In the experiment reported in Ref. the fused silica 
sample was exposed to a few-cycle pulse with a central 
photon energy ?icjl ~ 1-7 eV, and charge displaced by 
the optical field was collected by two gold electrodes (see 
Fig.II]). 
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Figure 1. A scheme of the metal-dielectric junction and the 
orientation of the laser field. Courtesy of A. Schiffrin. 



The dielectric slab is much thinner than the beam size 
100 /im) but much thicker than electron's de Broglie 
wavelength. In this case, effects of quantum confine- 
ment are negligible, and thus the stationary states are 
described by Bloch wavefunctions. An external field po- 
larizes the sample, and the induced surface charge creates 
a screening field. From the condition of continuity of the 
normal component of the electric displacement field at 
the surface between two media Di.^ = T)2.± it follows 
that the electric field E(t) inside the sample is connected 
with the laser field EL(i) and macroscopic polarization 
P(i) by the following relation: 



E(f) = EL(<)-4^P(i). 



(11) 



From the definition A = — E, we obtain the equation for 
the vector potential inside the dielectric: 



-A(t) = -EL(i)+47rP(t). (12) 
The macroscopic polarization satisfies the equation 
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where the charge current density is given by the following 
expressions^: 
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where averaging is done over a volume V, s and s' are 
spin indices of initial and final states, respectively, and 
Too is the free-electron mass. 
From this equation we obtain 



2e f d^k fv^ „ , 

BZ J'* 



X [5/,nk + Rep/,(k)] -iVvBeA(t)|, (15) 

where indices / and i enumerate all valence and conduc- 
tion bands, the sign function <Tfi is defined by Eq. ([SJ, 
the multiplier 2 accounts for the spin degeneracy, and 
Nyb is the number of valence bands. 

We solve Eqs. ([TOl), (US]), and ((13]) self-consistently for 
given laser field EL(i), energy spectrum ei(k), and mo- 
mentum matrix elements p/i(k). 

As long as the external field is not strong enough to 
generate a significant number of electron-hole pairs, the 
induced screening field can be evaluated from the linear 
polarization response: E -I- 47rP = eD. Regarding the 
response to be approximately instantaneous, we evaluate 
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the relative permittivity e at a central laser frequency wl 
and, from Eq. dTTl) . obtain 

A(i) - AlW/e. (16) 

In this case only Eqs. (fTO)) for the density matrix must 
be solved. 

In the next sections, we discuss the comparison of these 
two approaches and refer to the model that includes 
dielectric screening self-consistently [Eqs. ^TU\\ . (IT^ . 
and (fT5|) ] as OBE/SCDS, and to the instantaneous 
screening model [Eqs. (ITUl) and ([H])] just as OBE. 



III. RESULTS AND DISCUSSION 

It is well known that nonparabolicity of electronic 
bands plays an important role in strong-field phenom- 
ena^, and their description requires an accurate solution 
of the stationary problem for the entire Brillouin zone. 
We calculate the stationary electronic levels e^^k and mo- 
mentum matrix elements p/i(k) using the plane- wave 
pseudopotential methodi^>2S. 

The eigenvalue problem for the crystal is written as 



2mn 



C/(r) 



a,k(r) = er!,k(/>«,k(r), 



— Un.k(r) exp(ikr) is the wave function 
Wn,k(r) is a Bloch amplitude, and U{r) 



where 0„,k(i") 
of an electron 
is a periodical lattice potential. Since Un.kCr) has the 
periodicity of the lattice, it can be expanded in a Fourier 
sum: 



•n,k(r) = ^Ci'^exp(iK„,r), 



(17) 



where K„i are the reciprocal-lattice vectors. 

Assuming that the laser pulse is linearly polarized, we 
may reduce the problem to one spatial dimension and 
solve the density matrix equations for the k-states that 
belong to a certain direction in the Brillouin zone. In 
this paper, we consider a one-dimensional lattice with 
the pseudopotential 



U (z) — ci{l — tanh^ C2z), 
(n — l)a||/2 ^ z ^ na||/2, n € Z, 



(18) 



where a|| is the lattice parameter in the direction parallel 
to the laser field polarization. 

We assume that our pseudopotential describes the 
[001] direction of a-quartz, for which a|| = c = 5.4 A. 
Fit parameters ci = —2.2 and C2 = 0.9 were chosen such 
that the energy gap is Eg ^ 9 eV and the effective mass 
of an electron in the first conduction band is tUc ~ 0.4mo 
[see Fig. [5Ja)]. These values are in good agreement with 
the results for the F-A direction in the Brillouin zone of 
a-quartz obtained with more complex and rigorous treat- 
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Figure 2. (Color online) (a) Band structure obtained via 
numerical solution of the eigenvalue problem for the one- 
dimensional periodic potential given by Eq. (|18p . (b) Di- 
agonal matrix elements of the momentum operator Pa = 
hk + pu{k) versus k for the topmost valence and lowest con- 
duction bands. 



It is noteworthing that diagonal momentum matrix el- 
ements between the full Bloch functions Pa = hk+pu^k) 
quickly approach zero at the boundaries of the Brillouin 
zone, i.e. electrons are slowing down in this region [see 
Fig.EIb)]. 

In order to produce results that could be compared 
with experimental data, we must specify a proper connec- 
tion between our one-dimensional calculations and phys- 
ical values defined for a three-dimensional crystal. We 
assume that the electromagnetic wave is polarized along 
the z-axis and the current density slowly changes within 
the xy-plane. Then the current density J{t) in a 3D- 
crystal is connected with its ID counterpart via the ef- 
fective cross-section 6 of the three-dimensional Brillouin 
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where jiD{kz,t) is the one-dimensional analogue of the 
expression in braces from Eq. (jlSp , n is the unit vector in 
the direction of the current. In a first approximation, the 
parameter 6 can be estimated directly from the proper- 
ties of the lattice for a given material. For a-quartz with 
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a hexagonal lattice and current direction along the [001] 
direction, we get 



6si02 



3^/3 /27r 



1.2 at.u., 



where a± = a = 4.9 A is the lattice parameter of the 
crystal plane perpendicular to the direction of the laser 
field polarization. 

On the other hand, the value of & can be found more 
precisely from the requirement that OBE solution should 
give a correct linear response in a low-intensity region. 
With a given linear susceptibility x''^^ ^nd polarization 
calculated from Eq. (fTU)) . we have 6si02 ~ 1-43 at.u., 
which is in good agreement with the above estimation. 

For numerical integration of Eqs. (ITU)) . ([T^. and ([T^. 
we discretized the k-space and applied the fourth-order 
Runge-Kutta method with a constant stepsize. At the 
highest field intensities considered in this paper, numer- 
ical convergence was achieved with 18 bands and 201 k- 
points. The requirement of a large number of bands is 
a consequence of the velocity gauge drawbacks discussed 
earlier in Sec. [TTl 

In subsequent calculations, we use the following ex- 
pression for the vector potential of the laser field: 



Ai^{t) = -Ao6'(rL - |i|) cos"^ — sin(wLt + ^Pce) 
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where wl is the central frequency of the laser pulse, 0{x) 
is the Heaviside step function, Aq = Eq/wl and Eq are 
the amplitudes of vector potential and electric field in the 
vacuum, respectively, (pcE is the carrier-envelope phase. 
Total pulse duration tl is related to FWHM of |A(t)p 
as TL = 7rTFWHM/[4arccos(2-i/^)] sa 1.9tfwhm- 

Let us start the discussion of our numerical results with 
the distribution of populations in the first conduction 
band at the moment of time when the electromagnetic 
field is over. Fig. [3] shows that the few-cycle pulse gener- 
ates an asymmetric band population distribution in the 
crystal momentum space, which causes the appearance 
of macroscopic current. 

We argue that these symmetry-breaking effects of the 
waveform with a well-defined CE phase emerge due to the 
interference of different multiphoton excitation pathways. 
This phenomenon is schematically depicted on Fig. 21 
where arrows of different colors describe different mul- 
tiphoton channels that interfere with each other. Note, 
that excitation pathways might have the same number 
of photons with different frequencies, as well as different 
numbers of photons. In this picture, physical observables 
become CEP-dependent if the pulse is short enough, i. e. 
when its spectral width allows for existence of multipho- 
ton channels with odd and even numbers of photons for 
the same k-point in the Brillouin zone. The interference 
of optical excitation pathways'^'^ might be constructive for 
k and destructive for — k. When the laser pulse is over. 



electron-phonon collisions become the dominant interac- 
tion, which quickly restores the symmetry of electronic 
population distribution and make the current disappear. 

Figs. [2Ia) and (b) show that the maximum of pop- 
ulation distribution is shifted from the center, and the 
electron excitations spread over the whole Brillouin zone 
when the field intensity is increased. This situation is 
strongly opposite to the case of a weak field, where most 
of the charge carriers are situated around the extremal 
points in the Brillouin zone. At a relatively low field am- 
plitude oiEo^l V/A, both OBE and OBE/SCDS mod- 
els give very similar results, as expected, but for stronger 
fields, the second one shows the overall increase of the 
population and even more asymmetric distribution. This 
means that the screening field becomes sufficiently strong 
to cause interband transitions. 

For a deeper understanding of these observations, let 
us consider the polarization response. Fig. [SJa) depicts 
a typical time-dependent polarization in response to a 
laser pulse with a peak intensity near the damage thresh- 
old (Eq ^ 2 V/A). The polarization response within the 
first half of the pulse is linear because the field is not 
strong enough to excite charge carriers to the conduction 
band, so the contribution to the polarization comes only 
from small displacements of electronic shells that form 
the valence bands of the crystal. Once the field reaches 
a certain strength, electron-hole pairs are created, and 
the polarization response shows quantum beats with a 
characteristic energy approximately equal to band gap. 

A comparison of the plots in Fig.[5ja) shows that, with 
the self-consistently evaluated screening field which re- 
mains after the laser pulse, quantum beats are signifi- 
cantly reduced. The suppression of quantum beats in a 
strong field was already pointed out in our previous ear- 
lier paper™. In addition to that, the interaction of elec- 
trons with the screening field induced by quantum beats 
results in a phenomenon analogous to stimulated emis- 
sion. The induced screening field E(i) oscillates out of 
phase with respect to the polarization P(t), it drives tran- 
sitions from the conduction band to the valence band, 
and thus, decreases the amplitude of quantum beats. 

The time dependence of polarization [Fig. EJb)] shows 
that the self-consistent screening model predicts low- 
frequency oscillations that persist after the laser pulse. 
These oscillations occur as a result of collective electron 
motion driven by the surface charge field. For the pa- 
rameters of our simulations, the frequency of these os- 
cillations is in the terahertz region (e. g. 50 THz for 
Eq = 2.5 V/A), which is close to the frequency of plas- 
monic oscillations, although, in our model, we account 
for electron-electron interaction only implicitly, via the 
macroscopic screening field. In principle, these oscilla- 
tions should produce the electromagnetic radiation in 
the corresponding frequency range. The phenomenon 
of THz emission from cold plasma oscillations in semi- 
conductors excited by femtosecond optical pulses is well 
known and has been reported in a number of recent pa- 
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Figure 3. (Color online) Distributions of the lowest conduction band population ric after the laser pulse (Al ~ 800 nm, 
FWHM = 4 fs, 4>cE = 0) for different field amplitudes, calculated (a) with QBE and (b) with OBE/SCDS. Populations for 
positive and negative crystal momenta are shown with the dashed and solid lines, respectively. 




Crystal momentum 



Figure 4. (Color online) A schematic representation of multi- 
photon channel interference at different points in the Brillouin 
zone. 



pers^— . Also, there is an important analogy with the 
optical rectification effec t*^ ''i'^^ . 

Since we do not consider relaxation phenomena (radia- 
tive or non-radiative), these plasma oscillations do not 
decay with time in our simulations. Their amplitude 
might even grow, if the screening field is large enough 



to induce interband transitions. Also, the dielectric po- 
larization model does not take into account the energy 
loss due to the emission of electromagnetic radiation, so 
that the dielectric acts like a resonator for the induced 
currents. Consequently, the results of our present simu- 
lations are valid only within a small interval of a few fem- 
toseconds after the laser pulse. A more accurate descrip- 
tion of the strong-field phenomena requires the complete 
solution of semiconductor Maxwell-Bloch equations that 
take into account electron-electron and electron-phonon 
interactions, as well as electromagnetic radiation from 
accelerated charge carriers. 

In the experiment reported by SchifFrin et al. in Ref. [^, 
the measurements yielded the total transferred charge Qp 
in a fused silica junction defined as 

Qp = A,sqp, (20) 

where 



oo 



qp = P{t oo) = 




— oo 



is the transferred charge density, and A^s is the effective 
cross section of the active volume, estimated as ^ 5 x 
10-12 m^. 
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Figure 5. (Color online) (a) Polarization versus time calcu- 
lated with QBE (top panel) and with OBE/SCDS (bottom 
panel), Eq ^ 2 V/A, Al = 800 nm, ipcE ~ 0. (b) Polarization 
versus time for QBE (dotted line) and OBE/SCDS (solid hue) 
at the end of the laser pulse. For the field intensities close to 
the damage threshold, the low- frequency oscillations in polar- 
ization response are observed. 



This simple definition is directly applicable if theory 
includes all necessary relaxation mechanisms, both ra- 
diative and non-radiative, and the integral of charge cur- 
rent density over time takes a finite value. In order to 
compare our numerical results with the experiment, we 
need to estimate the value of integral (j2ip . taking into 
account the applicability limitations due to the absence 
of relaxation mechanisms in our OBE/SCDS model. As- 



suming that the current quickly decays after the laser 
pulse, we estimate the integral (j21l) by the value of po- 
larization right after the pulse. Thus, we can define the 
transferred charge density as an average value of the po- 
larization taken in the time interval larger than the pe- 
riod of quantum beats Tb and smaller than the period of 
the long-wave polarization oscillations Tp after the laser 
pulse at the highest considered field intensity 



^ J p{t)dt, Tb<Ai«rp 



In the subsequent calculations, we assume that At — 1 fs. 

In order to convert the calculated values from atomic 
to SI units, we use the following relation: 

qp [C/m^] = —jTT-^ X qp [at. u.] « 2.077 x qp [at. u.], 

where ae is the Bohr radius. 

Fig. [5] shows typical dependencies of the transferred 
charge density on the absolute value of carrier-envelope 
phase. The most striking difference between the results 
obtained with OBE (left panel) and OBE/SCDS (right 
panel) is the large shift (about OAtt) of the CE phase 
that maximizes the transferred charge. This phase shift 
has a simple explanation: in the model with constant 
screening the phase of the field inside the crystal is ex- 
actly the same as of the laser field. In the self-consistent 
screening model, the field inside the medium receives an 
additional phase shift from the polarization generated by 
the macroscopic charge current and quantum beats. 

As a justification of our model, we provide a compar- 
ison of our results with the measurements published in 
Ref. 0. Since a self-consistent evaluation of the screen- 
ing field is essential for very strong fields, we discuss only 
OBE/SCDS model. Fig.[7{a) shows very good agreement 
between theory and experiment up to Eq = 2.2 V/A. 
The discrepancy at high fields probably appears because 
the independent-particle approximation and dielectric 
screening model lose their applicability. 

At low field intensities, the major contribution to the 
transferred charge should come from the interference be- 
tween 5- and 6-photon channels. In this case, Qp is ex- 
pected to scale as E^'^^ , and we indeed find that the func- 
tion cx E^^ gives a good fit for Eq < 1.5 V/A [Fig. [T^b)]. 
For more intense pulses perturbation theory breaks down, 
and the measured transferred charge cannot be fitted by 
a simple power function. 

In order to verify the interpretation of the CEP- 
dependent currents as a consequence of the interference of 
multiphoton channels, we investigate the charge density 
as a function of pulse duration for fixed central frequency 
and peak intensity of the laser pulse [Fig.|5]. In the lower 
panel of Fig. [51 we depict the energies accessible via dif- 
ferent multiphoton channels. The shaded areas in this 
figure show the regions where the probability distribu- 
tion related to n-photon absorption 
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Figure 6. (Color online) Transferred charge density qp ver- 
sus carrier-envelope phase obtained from (a) OBE and (b) 
OBE/SCDS calculations with the following values of field am- 
plitude: 1.15, 1.2, 1.25, 1.3, and 1.35 V/A. 
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Figure 7. (Color online) (a) Intensity scan of the transferred 
charge maximized over CEP: comparison of OBE/SDCS cal- 
culations with experimental data, (b) Results of fitting the 
calculated transferred charge dependence (solid line) and ex- 
perimental data on a double-logarithmic plot with the power 
function of field amplitude. 
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Figure 8. (Color online) Top panel: CEP-optimized trans- 
ferred charge density versus FWHM of the laser pulse for high 
and low field amplitudes. Bottom panel: spectral overlap for 
5-, 6-, and 7-photon excitation probabilities as function of 
FWHM. 



exceeds 90% of its peak value. Here I{t) is the cycle- 
averaged intensity. 

Even though a longer pulse creates more charge carri- 
ers, the transferred charge decreases with the pulse du- 
ration for both low and high field amplitudes (Figs, 
and (b), respectively). This can be explained by the de- 
creasing overlap of multiphoton channels (see the bottom 
panel in Fig. |5]). 
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Figure 9. (Color online) Double-logarithmic plot of CEP- 
optimized transferred charge density and its fit with powers 
of Eo versus field amplitude for different central frequencies 
of the laser pulse (FWHM = 4 fs). 



Second, we consider the field amplitude scans of the 
transferred charge density for the pulses with different 
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central frequencies. Fig. IH] shows that numerically cal- 
culated field amplitude for the scans in the region 1- 
1.5 V/A may be fitted with power functions that 
correspond to the interference between n and n + 1 pho- 
ton channels. Despite the success of this interpretation, 
it is clear that this simple power law cannot be valid 
under all circumstances. For example, interference be- 
tween more than two multiphoton channels, transitions 
to higher conduction bands, and non-perturbative phe- 
nomena affect the scaling law. 

IV. CONCLUSIONS 

We have presented a quantum-kinetic approach for the 
description of strong-field non-resonant injection of pho- 
tocurrents in dielectrics and wide bandgap semiconduc- 
tors. In this model, the appearance of a non-zero trans- 
ferred electric charge and its dependence on the carrier- 
envelope phase are intrinsically related to an asymme- 
try of electronic population distribution in k-space. This 
asymmetry, which also remains after the laser pulse, can 
be explained by the interference of different multiphoton 
excitation channels. We have presented two arguments 
supporting this interpretation. First, we varied the band- 
width of the laser pulse and found that the amount of 
transferred charge is largely determined by the overlap 
of the dominant excitations channels. Second, we have 
shown that, for fields up to Eo ^ 1.5 V/A, the depen- 
dence of the transferred charge on the amplitude of the 



laser pulse is well described by the simple perturbative 
expression for interfering multiphoton excitation path- 
ways. The scaling law predicted by our numerical model 
is in very good agreement with the experimental data 
published in Ref. [a, except for the largest values of pulse 
amplitude Eq , where the concentration of excited charge 
carriers reaches a level where our mean-field description 
becomes inappropriate and plasma oscillations start to 
determine the polarization response. 

We have also compared two models for the evaluation 
of the screening field due to charges appearing on sur- 
faces of a mesoscopic structure exposed to a laser pulse: 
a model assuming an instantaneous linear dielectric re- 
sponse, and a more rigorous one where the polarization 
response is evaluated self-consistently with quantum dy- 
namics. We have found that the self-consistent evalua- 
tion of the screening field is essential for an accurate de- 
scription of CEP effects. For instance, these two models 
predict different dependencies of the transferred charge 
on the carrier-envelope phase and, thus, the different val- 
ues of (/)cE which maximize the current. 
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